Electronic and optical properties of ferromagnetic GaMnAs in a 
multi-band tight-binding approach 

Marko Turek/'Q Jens Siewert/ and Jaroslav Fabian^ 

^ Institut fiir Theoretische Physik, Universitdt Regensburg, D-93040 Regensburg, Germany 

(Dated: May 28, 2008) 

Abstract 

We consider the electronic properties of ferromagnetic bulk Gai_^Mn^As at zero temperature 
using two realistic tight-binding models, one due to Tang and Flatte and one due to Masek. In 
particular, we study the density of states, the Fermi energy, the inverse participation ratio, and 
the optical conductivity with varying impurity concentration x = 0.01 - 0.15. The results are very 
sensitive to the assumptions made for the on-site and hopping matrix elements of the Mn impurities. 
For low concentrations, x < 0.02, Masek's model shows only small deviations from the case of 
p-doped GaAs with increased number of holes while within Tang and Flatte's model an impurity- 
band forms. For higher concentrations x, Masek's model shows minor quantitative changes in the 
properties we studied while the results of the Tang and Flatte model exhibit qualitative changes 
including strong localization of eigenstates with energies close to the band edge. These differences 
between the two approaches are in particular visible in the optical conductivity, where Masek's 
model shows a Drude peak at zero frequency while no such peak is observed in Tang and Flatte's 
model. Interestingly, although the two models differ qualitatively the calculated effective optical 
masses of both models are similar within the range of 0.4 - 1.0 of the free electron mass. 

PACS numbers: 75.50.Pp, 71.55.Eq, 72.80.Ey, 75.50.Dd 
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I. INTRODUCTION 



Dilute magnetic semiconductors are investigated very actively for their potential of 
combining ferromagnetic with semiconducting properties. A prominent prototype is 
Gai_^Mn^As with Mn concentrations being typically x = 0.01 - O.IS.-*^'^'^ For these in- 
termediate to large concentrations critical temperatures up to 170 K could be observed.- 
The host material, GaAs, is a III-V semiconductor with a band gap of 1.5 eV. When the 
Mn impurities substitute the Ga atoms they act as acceptors which carry a local magnetic 
moment caused by the half-filled d-shell. This leads to a hole-mediated ferromagnetism.- 
The impurity binding energy is 0.11 eV.^ If the Mn concentration exceeds a critical value 
X ~ 0.01 the impurity wavefunctions at the Fermi energy overlap sufficiently for the material 
to undergo a transition towards a metallic state. 

A widely discussed question is whether the holes reside in an impurity band which is 
detached from the host valence band or in the valence band itself, which would cause dif- 
ferent transport properties.- To this end a variety of absorption experiments^'^'^'^ and 
measurements of the band gap and chemical potentiaU^i^ in GaMnAs were performed. For 
very low concentrations x <^ 0.01, an impurity band is formed with the Fermi energy re- 
siding therein.— This picture has also been used for the intermediate concentration range, 
X ~ 0.01 - 0.15, to explain optical measurements.^ On the other hand, there is experimental 
indication that the impurity band and the valence band have completely merged, see Ref. 
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and references therein. Another issue concerns the localization properties of the carriers. 
It has been argued that instead of the the impurity band interpretation the assumption of 
a merged valence and impurity band together with the existence of localized states in the 
band tail can explain the experiments.— 

Besides the experimental efforts for a better understanding of the electronic properties 
of Gai_^Mna.As a wide range of theoretical models has been developed.- Among these are 
first-principles calculationslS.tlS.^ effective single particle tight-binding approaches^*^*^"^, 
tight-binding approaches in combination with percolation theory2ii2^i2^^ dynamical mean 
field theories^i, effective theories based on k-p Hamiltonians2^i2Si'2Ii2^^ and large-scale Monte- 
Carlo studies of a real space Hamiltonian^^. The microscopic tight-binding approach has the 
advantage that one can study disordered systems in a non-perturbative way. The basis for 
this approach are the 16 sp'^ valence and conduction bands of GaAs that are approximated 
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very closely to experimental data throughout the entire Brillouin zone. System sizes up to 
approximately 2000 atoms can be treated within this approach with reasonable numerical 
effort. 

In this paper we present our results concerning the electronic structure of bulk 
Gai_^Mnj.As at zero temperature. We assume the system to be ferromagnetic by align- 
ing all Mn core spins into the z-direction. Specifically, we focus on the role of substitutional 
Mn impurities by neglecting interstitial disorder. The basis for our studies are well-known 
sp'^ tight-binding models for the host material GaAs which include spin-orbit coupling.— 
The disorder effects are treated in a non-perturbative way by changing the on-site and hop- 
ping terms of those Ga sites that are replaced by the Mn impurities. Instead of describing 
the effects of the Mn impurities by the two sp-d exchange constants A^^oa and Nq(3 onlyi^ii^ 
we follow two more elaborate approaches. The first approach, by Masek— (referred to as 
model A), uses first principles methods to determine the tight-binding parameters of the 
Mn impurities in a concentration regime around 10%. The second approach, by Tang and 
Flatte^"^ (model B), is based on a fit of the physically relevant tight-binding parameters to 
reproduce the binding energy of a single Mn impurity in the host material. Both approaches 
result in effective single particle tight-binding models in which the carrier-carrier correla- 
tions are included in a mean field way. This approximation is justified in the intermediate 
to high doping regime, i.e. x ^ 0.01.— Nevertheless, explicit inclusion of carrier-carrier 
interactions could lead to quantitative corrections which enhance localization effects.— 

In section [III we summarize the two approaches of how to find the reliable parameter sets 
for the inclusion of the Mn impurities and introduce our numerical method for solving the 
resulting equations in more detail. The explicit calculation of the eigenenergies and -states 
allows us to investigate the density of states and the position of the Fermi energy. In section 
mil we show that for Masek's model (A) the valence band is only slightly deformed while for 
the model by Tang and Platte (B) an impurity band forms which starts to merge with the 
host valence band at x ^ 0.01. The spatial extension of the states around the Fermi energy 
is investigated in section IIVI by means of the inverse participation ratio. We show that the 
states related to the impurity band show strong localization when their energy is close to the 
band edge. Based on the spectrum and the eigenfunctions we calculate the experimentally 
relevant optical conductivity. The relation of these results to the experimental findings is 
discussed in section |Vl 
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II. TIGHT-BINDING APPROACH 



Our analysis within the framework of a phenomenological tight-binding approacb^^i^ is 
based on an effective single particle Hamiltonian allowing for material-specific simulations 
of bulk GaMnAs systems. We considered disordered super-cells of up to 2000 atoms with 
the super-cell size being limited by the available computer resources only. In this work we 
focus on substitutional disorder by changing the on-site and hopping terms of certain Ga 
sites which are replaced by Mn. 

In order to introduce the terminology we briefly summarize the tight-binding approach. 
Using the single particle basis, 

)Mr-ta-R), (1) 
the eigenstates of the Hamiltonian can be written as 

o,a 

Here, o labels the atomic orbital 0o, cl specifies the atom sitting at the position ta within 
the super-cell, S is the total number of super-cells and R points to the several super-cells 
that are included. The sum over o and a in Eq. [2] thus runs over all orbitals of all atoms 
in a super-cell. Due to the periodicity and the Bloch character of the basis states Xkoa 
Hamiltonian is block-diagonal with respect to the different k vectors. The matrix elements 
in this basis thus read 

'H{k) 

where the sum goes over all the super-cells that contain the nearest or next to nearest 
neighbors, depending on the approximation employed. The vector f^n points from atom a 
to atom a' so that fVm = Ai?„„ + 1^/ — t^. On-site terms are thus characterized by fnn = 
while hopping terms correspond to f^m 0. Numerical diagonalization of the matrix ([3]) 
returns the eigenenergies and -vectors as a linear combination of the localized atomic orbitals 
(poa{r)- 

We have studied two different parameter sets for the inclusion of Mn impurities into 
existing tight-binding models for GaAs. The first set due to Masek (model A) is derived from 
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a first principles approach while the second one due to Tang and Flatte (model B) is based 
on phenomenologically deduced parameters. In either case we checked our simulation by first 
calculating the band structure and density of states for clean GaAs using the tetrahedron 
method.^- Then we considered disordered super-cells which are repeated periodically, and 
their ensembles. The sizes of the super-cells are typically between 128 and 1024 atomic sites 
while the disorder averages were performed over 5 to 15 random configurations, depending 
on the system size and the number of impurities. The finiteness of the systems limits 
the impurity concentrations to x ^ 0.006. As we will demonstrate below, the density of 
states is self-averaging for super-cells of the above-mentioned sizes so that increasing the 
number of disorder averages decreases the fluctuations of the results but does not give any 
new qualitative features. The self-averaging property also allows to imply periodic boundary 
conditions to the super-cells by which less fluctuating results can be obtained. These periodic 
repetitions of the super-cells correspond to a summation over various k points in the Brillouin 
zone. In our calculations we used between 20 and 505 different points in k space which are 
chosen in close resemblance to the ones given by the tetrahedron method. However, in 
contrast to the clean system, a linear interpolation of the band structure does not give an 
improvement of the quality of the results as the eigenenergies of a disordered super-cell for 
a given k lie too close to each other. 

1. Model A (Masek) 

The first model we study was introduced by Masek.— It starts from a tight-binding 
approach for the host material which includes the first and second nearest neighbor 
interactions'^^ as well as spin orbit coupling.— This gives a rather accurate description 
of the band structure including the conduction bands away from the F point. The on-site 
and hopping terms for the Mn impurity in the disordered system were obtained from a first 
principles calculation using a self-consistent Hartree-Fock-approximation for a single Mn 
impurity hybridized in a GaAs lattice. The resulting parameters were then checked within 
the same theoretical framework for a system with 10% of the Ga atoms replaced by Mn 
and appeared to be robust. This procedure gives for the on-site energies -Emh.sT = —0.4 eV, 
-^Mn.sj, = 0.0 eV, and Eu_n,-p = 4.374 eV. Additionally, the ten d-orbitals of Mn are exphc- 
itly included. These orbitals can be divided into two subgroups according to their sym- 
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metry, i.e., (Mn, t2) and (Mn, e). The corresponding values are i?Mn,t2T = —2.21981 eV, 
^Mn,t2i = 2.63427 eV, ^Mn.et = -3.01348 eV, and ^Mn.ei = 2.36445 eV. The hybridization 
of these d-levels with the s- and p-orbitals of the nearest neighbors is also included. The rele- 
vant hopping parameters are Vsd,a = -1-1077 eV, V^d,cr = —1.0341 eV, and I^^.tt = 0.4767 eV. 
This model is supported by photoemission spectroscopy experiments^- that revealed minor 
differences in the band structure of Gai-^Mn^^As compared to GaAs except for the appear- 
ance of the additional Mn d-states. 

2. Model B (Tang, Flatte) 

The second model we studied was introduced by Tang and Flatte.—*^ The electronic 
structure of the GaAs host material is described by a slightly different sp^ tight-binding 
model which also includes spin-orbit coupling but remains within the nearest neighbor 
approximation.— Despite the lack of the second nearest neighbor terms this model gives 
a fairly accurate description of the band structure of GaAs around the F point. The replace- 
ment of a Ga atom by a Mn impurity is modelled by an effective potential which describes 
the changes in on-site energies and hybridization of the Mn d-orbitals with the As p-orbitals. 
The Mn d-orbitals are included indirectly by a spin-dependent effective energy shift of the 
nearest neighbor p-orbitals. There are two independent parameters of this effective poten- 
tial. They are given by compensating the difference in the atomic ionization energies and 
by fitting the experimentally^ obtained Mn acceptor level lying at 0.11 eV. This model was 
applied to describe single Mn atoms and Mn pairs in GaAs and showed good agreement 
with experimental results for the local and total density of states as well as the shape of the 
wave functions.— 

III. DENSITY OF STATES 

The density of states, d{E, x, N), is obtained by direct diagonalization of the tight-binding 
Hamiltonian ([3]) and the evaluation of the sum 
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FIG. 1: Model A - Density of states for clean GaAs (black dotted, tetrahedron method) and 
disordered GaMnAs. The Mn concentration is x = 0.05. System sizes are N = 128 sites (505 k 
vectors, 15 disorder configurations) and N = 1024 sites (20 k vectors, 5 disorder configm'ations). 
The arrow indicates the position of the Fermi energy, Ep = —0.34 eV. The upper (lower) part of 
the figure shows the spin up (down) contribution to the density of states. 
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FIG. 2: Model B - Density of states for clean GaAs (black dotted, tetrahedron method) and 
disordered GaMnAs. The Mn concentration is set to x = 0.05. System sizes are = 128 sites (505 
k vectors, 15 disorder configurations) and N = 1024 sites (20 k vectors, 5 disorder configurations). 
The arrow indicates the position of the Fermi energy, Ep = 0.30 eV. The upper (lower) part of the 
figure shows the spin up (down) contribution to the density of states. 



where Se{E) is a broadened delta-function of the chosen width e = 5 meV. The eigenvalues 
E^^ of the Hamiltonian are labeled by the band index b. The density of states is normalized 
by the total number of states n{oo). As the parameter sets for clean GaAs are different for 
the two models the resulting densities of states also show different behavior, in particular 
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FIG. 3: Model A - Density of states of GaMnAs. The Mn concentration is set to x = 0.009 
(1024 atoms, 20 k vectors, 5 disorder configurations), x = 0.051 (432 atoms, 89 k vectors, 10 
disorder configurations), and x = 0.148 (432 atoms, 89 k vectors, 10 disorder configurations). The 
arrows indicate the positions of the corresponding Fermi energies: Ep{x = 0.009) = —0.17 eV, 
Ep{x = 0.051) = —0.34 eV, E'f{x = 0.148) = —0.49 eV. For comparison the densities of states for 
clean GaAs and MnAs are shown. 

for the conduction band, see Figs. [1] and [2l 

The results for a Mn concentration x = 0.05 and two different system sizes are shown in 
Figs. [Hand [2] for the models A (Masek) and B (Tang, Flatte), respectively. In both cases 
one can clearly observe the self-averaging properties of the density of states - the deviations 
between the two systems of different sizes are on the same order as the fluctuations due to 
the disorder. This self-averaging property of the density of states can be observed for the 
entire range of concentrations studied. This implies that the quantity d{E, x, N) = d{E, x) 
is independent of the system size if the smaller fluctuations are neglected. 

For model A (Masek) one finds that the density of states is only slightly perturbed in 
comparison to the density of states for the clean system. This is in contrast to the results 
of model B (Tang, Flatte) where, especially around the top of the valence band of the 
host material, the number of spin-up states is significantly enhanced. In either model the 
spin-down contribution of the valence band remains unaltered. 

As energies around the Fermi energy are most interesting we present the density of states 
for a smaller energy window and various Mn concentrations in the Figs.[3]and|H For model A 
(Masek), increasing the Mn concentration from 1% to 15% does not change the qualitative 
but only quantitative characteristics of the density of states. In particular, no detached 
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FIG. 4: Model B - Density of states of GaMnAs. The Mn concentration is set to x = 0.009 (686 
atoms), X = 0.051 (432 atoms), and x = 0.148 (432 atoms). In either case 89 k vectors and 10 
disorder configurations were used. The arrows indicate the positions of the corresponding Fermi 
energies: (0.009) = 0.10 eV, (0.051) = 0.30 eV, (0.148) = 0.65 eV. For comparison the 
densities of states for clean GaAs and MnAs are shown. 

impurity band forms within this model and the given range of Mn concentrations. The size 
of the band gap decreases as some of the spin-up states are shifted from below the top of 
the valence band into the gap of the host material. However, as adding Mn impurities also 
adds holes, the Fermi energy decreases and moves deeper into the valence band, see also 
Fig.O Thus, the major effect of an increasing number of Mn impurities is the larger number 
of holes which leads to a lowering of the Fermi energy. This trend qualitatively follows the 
Fermi energy of pure GaAs with the corresponding number of holes added, see Fig. |6l 

For model B (Tang, Flatte) a completely different picture arises. Due to the strong 
effective potential of the Mn impurities there is a large number of states shifted from the 
host valence band far into the gap. Already at concentrations of 5% there is almost no gap 
left at all. The second effect of the addition of Mn impurities is again the increased number 
of holes which lowers the Fermi energy. However, for model B (Tang, Flatte) this increasing 
number of holes does not compensate the effect of the large number of states appearing 
in the host gap. Hence one observes an increasing Fermi energy, measured relative to the 
top of the host valence band, see Fig. [51 The second major difference to model A (Masek) 
is the appearance of a detached impurity band for small enough impurity concentrations. 
This impurity band starts to merge with the valence band for concentrations around 1%, see 
Fig. m which is consistent with experimental data.-^S Its center is positioned at the impurity 
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FIG. 5: Model A, B - Fermi energy of GaMnAs as a function of the impurity concentration 
measured from the top of the GaAs host valence band. The Fermi energy of model A moves into 
the host valence band while it increases into the host gap for model B. System sizes between 128 
and 1024 atoms are shown. The dotted lines serve as a guide to the eye. 
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FIG. 6: Model A, B - Fermi energy measured relative to the top of the GaMnAs valence band 
and gap size Agap. The dotted lines serve as a guide to the eye. The dashed line shows the Fermi 
energy of p-GaAs without impurities where the number of holes is adjusted to the corresponding 
GaMnAs impurity concentration x. 

binding energy of Mn while its half width is on the order of 0.1 eV. For x ^ 0.05 one cannot 
distinguish the impurity band from the valence band anymore as the two bands have merged 
completely. 

The results for the Fermi energy and the band gap are summarized in Figs. [5] and [61 
In Fig. [5] the Fermi energy of GaMnAs measured with respect to the top of the GaAs 
host valence band is shown. From the behavior of the Fermi energy as a function of the 
impurity concentration it is evident that the two effects of adding holes and shifting states 
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are of unequal importance for the two models, as described in the previous paragraphs. In 
Fig. [6] we present the Fermi energy measured with respect to the GaMnAs valence band. By 
choosing these concentration dependent reference points for the energy we partially suppress 
the effects due to the shifted states. Hence, the Fermi energy decreases with respect to this 
energy reference point for both models. Again one finds that the addition of the holes is the 
major effect for model A (Masek) as the resulting Fermi energy behaves very similar to the 
case of clean p-GaAs with the corresponding number of holes added, see Fig. |6l Furthermore 
we have included the results for the band gap Agap into Fig. [61 This information can be used 
in the design and study of spintronic devices including hole injection into GaMnAs. The 
stronger influence of the disorder in model B (Tang, Flatte) compared to model A (Masek) 
becomes evident as the size of the gap is much faster decreasing for model B. Experiments 
based on scanning tunneling microscopy conductance spectra resulted in an estimate Agap ~ 
1.23 eV for x ~ 0.03 (Ref. Il2l ) and, more recently, photoconductivity measurements on 
GaMnAs heterostructures yielded Agap ~ 1.41 eV for the same Mn concentration.— The 
chemical potential was found to be in the vincinity of the Mn impuritiy level.— However, the 
considered samples were unannealed or annealed only for short times implying a significantly 
reduced carrier concentration due to Mn intersticials. Therefore, a direct comparison with 
our data is not possible and further, more detailed, experimental investigations are necessary 
in order to decide which model provides a more accurate description. 

The differences for the density of states between the two models can be better understood 
when the density for a clean MnAs (x = 1) system is calculated. As zinc-blend MnAs does 
not form a stable configuration^ this calculation does not correspond to a real physical 
system. Nevertheless it describes the x ^ 1 limit of the two models giving a good indication 
of what to expect for increasing x. The parameterization of model A (Masek) results in an 
increased number of states for MnAs in the gap of GaAs. For the disordered GaMnAs there 
is a clear tendency for the density of states to change from the GaAs-shape (x = 0) towards 
the MnAs shape (x = 1). This can, for example, be seen by the formation of the shoulder in 
the density of states &t E —\ eV and x = 0.148, see Fig. [3l For model B (Tang, Flatte), 
the MnAs density shows a qualitatively new feature a.t E ^ 0.5 eV in the gap of the GaAs, 
seen in Fig. |H The appearing peak is due to the shift in the energies of the p-states which 
is reflected in the rapidly vanishing gap for the disordered GaMnAs. 
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FIG. 7: Model A - Scaled inverse participation ratio [N-I{E, x, N)] for fixed concentration x = 0.05 
and various system sizes N = 128 (505 k vectors, 15 configurations), N = 432 (89 k vectors, 10 
configurations), and N = 1024 (20 k vectors, 5 configurations). The arrow indicates the position 
of the Fermi energy Ep = —0.34 eV. 

IV. INVERSE PARTICIPATION RATIO 



Besides the density of states and the Fermi energy we also studied the character of the 
eigenstates, in particular their localization properties. To this end we analyzed the inverse 
participation ratio averaged over all states in a given energy window: 
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HE, X, N) = '-^ y . (5) 
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b,k 



Qualitatively, the participation ratio corresponds to the average number A^^occ of sites occu- 
pied by the eigenstates of given energy E, i.e. 1// ~ N^cc- Thus, I{E,x,N) scales with 
the inverse system size for extended states as the probability Ic^oaP to find the electron 
at a certain atom is on the order of the inverse system size. For strongly localized states 
I{E,x,N) remains constant independently of the system size. 

First, we consider the scaling of [A^ ■ I{E, x, N)] with system size for a fixed impurity 
concentration. The results for the two models are shown in Figs. [7| and M where [A^ • 
I{E,x = 0.05, A^)] is plotted. As expected one finds extended states for the energies deep 
within the valence band. As the energy is increased towards the Fermi energy, the value of 
[N-I{E, X, N)] increases, meaning that on average less atoms participate in the corresponding 
states. Comparing Figs. [7] and [8] one finds this effect to be much stronger for model B (Tang, 
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FIG. 8: Model B - Scaled inverse participation ratio [N-I{E, x, N)] for fixed concentration x = 0.05 
and various system sizes N = 128 (505 k vectors, 15 configurations), N = 432 (89 k vectors, 10 
configurations), and N = 1024 (20 k vectors, 5 configurations). The arrow indicates the position 
of the Fermi energy Ep = 0.30 eV. 

Platte). Nevertheless a clear linear scaling of [A^ ■ I{E-f,x,N)] with system size N, which 
would indicate strongly localized states at the Fermi energy, cannot be found, see Fig. [9l 
This implies that the states around the Fermi energy reside on a subset of lattice atoms; yet 
they are still extended throughout the sample. This holds true for concentrations between 
5% and 15%, see Fig. [91 For model B (Tang, Flatte), this effect is more pronounced and we 
will argue that the states are spread over the sub-lattice of Mn sites. As the density of states 
decreases for energies closer to the band edge there are less states giving a contribution to 
the inverse participation ratio ([5]). This causes the stronger fluctuations in I{E,x,N) that 
can be observed for energies around and above the Fermi energy. 

Going from the Fermi energy Ep to higher energies close to the band edge, -Ebe? the states 
remain delocalized for model A (Masek) while they show strong localization for model B 
(Tang, Flatte). Qualitatively this conclusion is supported by the following argument. If an 
impurity state is strongly localized it cannot overlap with other impurity states located at 
neighboring impurity sites. This implies that the number of sites occupied by an impurity 
state must be significantly smaller than N/Nun- As the participation ratio is a measure for 
the number of occupied sites N^cc it follows that the states are strongly localized only if 
I{E,x,N) is large enough so that x ^ I{E,x,N) independent of the system size A^. This 
is consistent with the statement that I{E,x,N) scales inversely with the system size for 
extended states while it remains constant for localized states. For model A (Masek) the 
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condition x ^ I{E,x,N) is never fullfilled in the entire energy range, see Fig. [71 However, 
for model B (Tang, Flatte), going beyond the Fermi energy this criterion can be fullfilled 
implying that there are localized states, see Fig. [HI 

To demonstrate this localization property for model B (Tang, Flatte) in more detail we 
show the maximum values of [A^ ■ I{E, x, N)], which correspond to energies i?be close to the 
band edge, in Fig. [TOl The clear linear increase of the scaled inverse participation ratio for 
the states closer to the band edge indicates strong localization at E around i?be- On the other 
hand, at lower energies E = Ep+0.5 eV a transition from delocalized to localized states takes 
place depending on the concentration of Mn impurities, see the inset of Fig. [TOl The linear 
increase of [A^ ■ I{Ep + 0.5 eV, x, N)] with size indicates localized states for concentrations 
X < 0.15. On the other hand, the scaled inverse participation ratio remains constant for x = 
0.15 for all system sizes meaning that the states are extended. Obviously, the concentration 
of Mn impurities is in this case large enough for the impurity states to overlap. Generally one 
can assume that for each fixed energy E there is a critical concentration Xcnt{E) for which 
the slope of [A^ ■ I{E, x, N)] vanishes. Based on this critical concentration one can estimate 
the spatial extension of the Mn impurity states. The average distance of impurities for given 
concentration x is approximately 0.4 nm ■ x~^^^. For the critical concentration Xcrit(-£') the 
impurity states start to overlap as the transition from localization to delocalization takes 
place. Therefore the spatial extension of the states at a fixed energy can be estimated by 
0.4 nm ■ (E). From Fig. [8] we find Xcrit(-£'F + 0.5eV) ~ 0.15 implying that the impurity 
states extend over a range of ~ 0.7 nm. At higher energies, i?be, extrapolation of our data 
leads to Xcrit(-Ebe) ~ 0.18 and a spatial extension of the impurity states of ~ 0.6 nm. This 
estimate is consistent with the experimentally deduced value.— 

In a second step we analyze the infiuence of the Mn concentration on the inverse partic- 
ipation ratio ([5]). Figures [TT] and [T2] show the results for a fixed system size A^ = 432 for 
models A (Masek) and B (Tang, Flatte), respectively. For comparison, we have also included 
the result of a larger system with a smaller concentration of Mn impurities. In the case of 
model A (Masek), the increase of [A^ ■ I{E, x, N)] with increasing energy starts first for the 
highest concentration x = 0.15. The reason can be found in the density of states - the larger 
X is, the stronger the valence band changes. However, if one considers [A^ ■ I{Ep,x,N)] at 
the Fermi energy, which also depends on x, the infiuence of the concentration is negligible, 
see Fig. [T31 where 1/[A^ ■ /(-Ef, N)] is plotted as a function of x. 
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FIG. 9: Model A, B - Scaled inverse participation ratio [N ■ I{E-p, x, N)] at the Fermi energy as a 
function of the system size A^. Mn concentrations x = 0.05, 0.10, and 0.15 are shown for model A 
and B. 
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FIG. 10: Model B - Maximum value of the scaled inverse participation ratio [N ■ I^i^^{E\^c,x,N)] 
for energies close to the valence band edge as a function of the system size N. The inset shows the 
same quantity at energies Ep + 0.5 eV. 

For model B (Tang, Flatte) the situation is different since [A^ ■ I{Ep, x, N)] increases with 
decreasing concentration x, as shown in Fig. [12], meaning that the number of occupied sites 
increases with x. This suggests that the states around the Fermi energy tend to the Mn 
sub-lattice. This interpretation is supported by the fact that 1/[A^ ■ I{Ef,x,N)] increases 
linearly with the concentration x, see Fig. [131 As 1/[A^ ■ I{Ef,x,N)] is a measure for the 
relative number of atoms where the states reside one can conclude that for model B (Tang, 
Flatte) the states tend to spread on the Mn sub-lattice and nearby sites. As an independent 
check of this conclusion we further looked at the probability of states with given energy to be 
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FIG. 11: Model A - Scaled inverse participation ratio \N ■ I{E,x, N)] for various concentrations 
X = 0.009 (1024 atoms, 20 k vectors, 5 disorder configurations), x = 0.051, 0.102, and 0.148 (432 
atoms, 89 k vectors, 10 disorder configurations). The arrows indicate the positions of the Fermi 
energy. 
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FIG. 12: Model B - Scaled inverse participation ratio \N ■ I{E,x,N)] for various concentrations 
X = 0.009 (1024 atoms, 20 k vectors, 5 disorder configurations), x = 0.051, 0.102, and 0.148 (432 
atoms, 89 k vectors, 10 disorder configurations). The arrows indicate the positions of the Fermi 
energy. 

found on the Mn impurities. As a result (not presented here) we found that this probability 
is indeed strongly increased already at the Fermi energy. It increases to almost one at even 
higher energies i?be close to the band edge. 
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FIG. 13: Model A, B - Scaled participation ratio 1/[N ■ I{Ep,x,N)] at the Fermi energy as a 
function of the Mn concentration x. System sizes of iV = 1024, 686, and 432 atoms are shown for 
model A and B. 

V. CONDUCTIVITY 

The optical conductivity is a directly measurable quantity which is closely related to the 
absorption. As there have been various experiments studying this quantity for GaMnAs we 
analyze the conductivity obtainable from the two models A (Masek) and B (Tang, Flatte) 
and compare our results with the experimental findings. The starting point for our analysis 
is the expression for the real part of the conductivity. 



m^A ^ J (27r)3 

BZ 



d^A; 



(/) 



(6) 



where / (i) labels the final (initial) states and pz is the 2;-component of the momentum 
operator. We restrict our considerations to a linearly polarized field in z-direction. As we 
investigate the optical conductivity in the infrared regime the eigenstates must share the 
same wavevector k. The resulting conductivity cannot be expected to be quantitatively exact 
as tight-binding approaches are known to underestimate the absorption.— Nevertheless, 
the order of magnitude of various related quantities, like effective masses, and important 
qualitative conclusions can be extracted. Within the tight-binding approach the matrix 
element of the momentum operator in Eq. E] can be expressed in terms of the Hamiltonian 
matrix, Eq. [31 and the distance between the localized orbitals.— This approximation neglects 
the spin-orbit coupling for the evaluation of the matrix element itself but correlations of the 
density of states are treated correctly giving a reasonable result for the optical conductivity.— 
The Figures [H] and [15] show the conductivity obtained from models A (Masek) and B 
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FIG. 14: Model A - Conductivity Re(T(A,x) for various concentrations x = 0.009 (686 atoms, 89 
k vectors) and x = 0.051, 0.101, and 0.148 (432 atoms, 240 k vectors). The results are averaged 
over 8 different disorder configurations. The inset shows the results for p-doped GaAs in the clean 
limit of the model with the corresponding number of holes added. The dashed-dotted line is a 
Lorentzian fit to the Drude peak at A = for the x = 0.051 data. 
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FIG. 15: Model B - Conductivity Recr(A,x) for various concentrations x = 0.009 (686 atoms, 89 
k vectors) and x = 0.051, 0.101, and 0.148 (432 atoms, 240 k vectors). The results are averaged 
over 8 different disorder configurations. The inset shows the results for p-doped GaAs in the clean 
limit of the model with the corresponding number of holes added. 

(Tang, Flatte), respectively. In either case, the low frequency values are of limited accuracy 
as the system is of a finite size which implies a finite minimal level spacing. The minimal 
cut-off energy difference is set here as Amin ~ 5 meV by using a broadened 5-function of 
this width in Eq. [6l 

As there is a clear qualitative difference between the two models let us start our discussion 
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with model A (Masek). For this model we find a low- frequency Drude peak which has 
also been observed in experiments.- For low concentrations, x ^ 0.01, additional peaks at 
~ 0.15 eV and ~ 0.33 eV can be identified clearly. These additional peaks correspond to 
inter- valence band transitions which blue-shift to higher energy differences and become wider 
as the number of impurities increases. This is very similar to the experimental observations 
in p-doped GaAs.-^ Eventually, for concentrations x ^ 0.05, these peaks are vanishing in 
the background of the Drude peak centered at A = 0. The width of this Drude peak can be 
related to a scattering time r. This width and the peak height can be estimated by fitting a 
Lorentzian to the conductivity curve for small A,- see Fig. [TH As this fitting procedure can 
only be performed with some uncertainty the resulting scattering times should be considered 
as qualitative estimates. For the concentrations studied they range from r ^ 27 fs to 6 fs, see 
Table [H The saturation of the scattering time at high concentration is due to the break down 
of the independent scattering mechanism. The height of the peak gives the DC-conductivity, 
e.g., (Tdc ~ 4000 Q^^cm^^ for x = 0.05, which is one order of magnitude larger than the 
experimentally observed values.-"^!^ Additionally the peak height and width are related to 
an effective mass m*. The estimation of this effective mass is based on a sum rule for the 
evaluation of the plasma frequency. 



with 77-0 being the carrier density and A^ax the cut-off energy. When Eq. [7] is applied to the 
Lorentzian fit of the Drude peak the cut-off energy can be set to infinity. This procedure 
of fitting a Lorentzian and applying the sum rule leads to effective masses m* /nif. ^ 0.9, 
see Table HIl Based on the estimates for the scattering times and these effective masses 
we conclude that {kpl)* ^ 10 - 14, see Table [H for x = 0.01 - 0.15 which means that the 
assumption of weak scattering is applicable. However, the compensation of some of the 
Mn acceptors due to intersticials lowers the effective hole density and thus also decreases 
the values for {kpl)* accordingly. The estimates based on experimental measurements lie 
between kpl ^ 0.75^^'^^ and kpl ^ 1 - 5^'^. Finally, we give an estimate of the mobility 
/i ^ 13 cm^/Vs for x = 0.1 which is by a factor 3 — 5 larger than the experimental values.-^ii 
A different approach based on the sum rule is to apply Eq. [7] directly to the numerical 
data and choose a cut-off that includes all inter-band transitions but excludes valence band 
to conduction band transitions^, i.e. A^ax = 1-5 eV. This means that not only the narrow 
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peak at A = is considered but the wider overall shape of cr(A) which decays on an energy 
scale ~ 0.15 eV. This defines an optical mass mopt which, for model A (Masek), remains 
approximately constant, i.e. mo^^/irtf. ^ 0.5, for x = 0.01 - 0.15, see Table [ITl This is to be 
compared with a stronger increase of the magnitude of the optical masses niopt/^e = 0.6 - 
1.2 which we obtain by the same method for p-doped GaAs with the corresponding number 
of holes added. This discrepancy not only in the magnitude but also in the rate of the 
increase of rriopt/me with respect to the concentration x indicates that the band structure 
of the host valence band is altered due to the disorder. 

In strong contrast to model A (Masek), there is no low- frequency Drude peak for model 
B (Tang, Flatte). For low concentrations x = 0.009 the impurity band causes an absorption 
maximum around 0.1 eV. Its half-width of ~ 50 meV refiects the width of the impurity 
band. For higher concentrations the impurity band completely merges with the host valence 
band leading to a more or less constant optical absorption. This behavior is in qualitative 
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At concentrations x ~ 0.1 a 



agreement with the experimental observations in Refs. |9| and 
broad feature with a maximum at 0.5 eV forms due to the larger number of impurity states 
which are shifted into the host gap. The absolute values of the conductivity lie between 
100 Q^^cm^^ and 900 Q^^cm^^ for the various concentrations which is in good agreement 
with the measurements.-' Applying the sum rule ([7j) with A^ax = 1.5 eV which defines the 
optical mass we find mopt/iTie ~ 0.4 for x = 0.009 increasing to iTLopt/rrie ^ 0.9 for x = 0.148, 
see Table HIl The corresponding values of p-doped GaAs range from mopt/w^e ~ 0.6 to 
fTiopt/^e ~ 1.2. As we have described in the previous section, model B (Tang, Flatte) is 
close to a metal-insulator transition for the impurity concentrations considered. This is 
refiected by the very broad peak implying very short scattering times which are inconsistent 
with the weak scattering limit. For the larger concentrations a; ^ 0.1 we estimate the half- 
width of the broad maximum of Re(T(A) to be ~ 0.5 eV implying r ~ 1 fs. From this 
follows that /i ~ 3 cm^/Vs and /cf / ~ 2.5 so that the mean free path is on the order of the 
wavelength of the carriers. Hence the loffe-Regel-limit is reached and a Drude-Boltzmann 
description of the conductivity cannot be expected to be a valid approximation any more. 

In Ref. isl the reasoning for the existence of an impurity band was based on the analysis 
of effective masses. However, the same analysis performed for our simulation data does not 
give a useful tool to distinguish the case without an impurity band from the one where an 
impurity band can still be identified. Although the density of states and the localization 
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Mn concentration x 


0.01 


0.05 


0.1 


0.15 


model A: rlfs] 


27 


7 


6 


6 


model A: (kpl)* 


10 


14 


14 


14 


model A: (fcpOopt 


6 


17 


18 


14 


model B: (AjpOopt 






2.5 


2.8 



TABLE I: Estimates for the dimensionless scattering length kpl. The values for (A;f/)* are based 
on a Lorentzian fit while (A;F/)opt was obtained by direct integration of the numerical data. The 
{kpl) values for model B indicate that for the given range of parameters the Drude-Boltzmann 
theory is not unambiguously applicable. We have nevertheless included the values for comparison 
with experimental data. 



Mn concentration x 


0.01 


0.05 


0.1 


0.15 


model A: m* /me 


1.0 


0.5 


0.9 


1.0 


model A: mopt/'m-e 


0.5 


0.4 


0.5 


0.5 


model B: mopt/m-e 


0.4 


0.6 


0.7 


0.9 



TABLE II: Estimates for the effective masses based on the sum rule ([7]). The values for m* /me 
have been obtained by the Lorentzian fit to the Drude peak obtained for model A (Masek). The 
masses Jriopt/w-e were calculated by integrating the numerical data directly. 

properties are very different for the two models A and B we find rather similar values for 
the effective and optical masses, see Table [Tll 

VI. CONCLUSIONS 

We have presented tight-binding studies of the electronic and optical properties of 
Gai_a;Mna,As based on two different parameterizations. The first parameterization due to 
Masek (model A), is based on first-principles calculations. In this model, the Mn impurity 
is described by a change in the on-site and hopping terms as well as the additional inclu- 
sion of the d-orbitals. This approach leads to certain qualitative changes in the electronic 
properties, i.e., the density of states. However, it does not give rise to the formation of 
an impurity band. Instead, the inclusion of the disorder has a very similar effect as the 
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inclusion of an equally large amount of holes without additional disorder. For the optical 
properties that can be deduced from Masek's model (A), we find the formation of a Drude 
peak as a clear signature of the disorder. This is in agreement with recent experiments. The 
estimations for the dimensionless scattering length {kp I) and the mobility fi give reasonable 
results while the absolute value of the DC conductivity is significantly larger than those 
observed in experiments. The inter-band transitions, that are also identifiable for low con- 
centrations x < 0.1, blue-shift with increasing amount of disorder which is in contradiction 
to the experiments performed on GaMnAs. 

The second parameter set, by Tang and Flatte (model B), is based on a more phe- 
nomenological modelling of the Mn impurities. It describes the bound level of a single Mn 
impurity rather well and leads to a formation of an impurity band in the density of states 
of Gai_2;Mna;As for small x ^ 0.01. This impurity band merges with the host valence band 
at concentrations x ~ 0.01. The Fermi energy lies within this impurity band and the states 
at the Fermi level are extended for Mn concentrations between 5% and 15%. However, the 
eigenstates become localized for higher energies closer to the band edge. The optical con- 
ductivity is characterized by the absence of a Drude peak at zero frequency and a rather 
featureless shape. It shows a broad maximum for larger Mn concentrations with a half-width 
of approximately 0.5 eV. From this width we deduce short scattering times leading to the 
conclusion that a Drude-Boltzmann theory of weak scattering is not applicable. 

As the experiments on the electronic properties of GaMnAs are not entirely conclusive 
at the moment it is difficult to judge which model fits best for a given physical quantity. 
Either model describes some of the experimental findings correctly while it fails for others. 
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